Damping of Bogoliubov Excitations in Optical Lattices 
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Extending recent work to finite temperatures, we calculate the Landau damping of a Bogoliubov 
excitation in an optical lattice, due to coupling to a thermal cloud of such excitations. For simplicity, 
we consider a ID Bose-Hubbard model and restrict ourselves to the first energy band. For energy 
conservation to be satisfied, the excitations in the collision processes must exhibit "anomalous 
dispersion" , analogous to phonons in superfluid 4 He. This leads to the disappearance of all damping 
processes when Un c0 > 6t, where U is the on-site interaction, t is the hopping matrix element and 
n c0 (T) is the number of condensate atoms at a lattice site. This phenomenon also occurs in 2D and 
3D optical lattices. The disappearance of Beliaev damping above a threshold wavevector is noted. 
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Recently, several experimental papers have reported 
results on the damping of collective modes in Bose con- 
densates in a one-dimensional (ID) periodic optical lat- 
tice potential Q. Most theoretical studies [E IE IE IE 
[E 0| have concentrated on the case of a pure conden- 
sate at T = and have ignored the effect of a ther- 
mal cloud, which can cause damping of condensate ex- 
citations at finite T. We discuss a Bose-Hubbard tight- 
binding model using the Gross-Pitaevskii (GP) approxi- 
mation, but generalized to include a static thermal cloud 
of non-condensate atoms. For illustration, we consider 
a ID model but emphasize the same features appear in 
2D and 3D optical lattices [3j. We limit our analysis to 
relatively strong optical lattices such that the lowest en- 
ergy band is fairly narrow. We calculate the temperature 
dependence of the number of condensate atoms n c0 (T) in 
each lattice well, as well as the Landau damping of con- 
densate modes due to coupling to thermal excitations. 
We find that for damping processes to occur, the dis- 
persion relation E q of the Bogoliubov Bloch excitations 
must initially bend upward as the quasi-momentum q in- 
creases. This "anomalous dispersion" is also the source 
of 3-phonon damping in superfluid 4 He 0, . This con- 
dition leads to the disappearance of all damping in an ID 
optical lattice when the normalized interaction strength 
a = Un c0 /t > 6 (where U is the on-site interaction and 
t is the hopping matrix element). More generally, all 
damping processes involving three excitations vanish in 
a /^-dimensional simple cubic tight-binding optical lat- 
tice when a > 6D [E . This ability to control the excita- 
tion damping in optical lattices may be very important 
in applications. 

Since our major interest is in the thermal gas of ex- 
citations in an optical lattice, it is important to keep in 
mind the distinction between the Bloch excitations as- 
sociated with linearized fluctuations of an equilibrium 
Bose-condensate and the stationary states of the time- 
independent GP equation for the Bose order parameter. 
In a continuum model, the latter states [E Q can be de- 



scribed by the eigenfunctions $®(x) = e lkx Uk(x), where 
the Bloch condensate function satisfies the usual peri- 
odicity condition u^ix) = Uk(x + Id), where d is the 
optical lattice spacing and I is an integer. Physically, 
3>2(2:) corresponds to a solution of the static GP equa- 
tion with a superfluid flow in the periodic potential, with 
the condensate quasi-momentum being given by hk. The 
fluctuations of these states are also described by a quasi- 
momentum q in the first Brillouin zone (BZ) and will 
be referred to as Bogoliubov Bloch excitations of the op- 
tical lattice. The thermal cloud at finite temperatures 
is an incoherent gas of these Bogoliubov excitations in 
the first band. In our tight-binding model, the Bloch 
condensate function is = e lkld ^J n c0 and we only 

consider a Bose condensate in the k = Bloch state, i.e., 

The ID optical lattice model (along the z axis) has 
been discussed in many recent papers at T = [E IE 
IE El- The optical lattice is described by V op (z) = 
sEr cos 2 (fcz), where Er = ^Mj- is the photon recoil en- 
ergy. We assume that the band gap energy between the 
first and the second excitation band is large compared 
to the temperature (2k^,T fE^ -C s), and thus only the 
first band is thermally occupied. The radial trapping 
frequency is assumed to be very large, so that motion 
in this direction is effectively frozen. The depth of the 
optical lattice wells is assumed large enough to make the 
atomic wave functions well localized on the individual 
lattice sites, described by a tight-binding approximation. 
This Hamiltonian is the Bose-Hubbard model, 

H = -t^(a, T +1 a, +a]a l+1 ) + -U^ajaja^, (1) 
i i 

where a/ and a, are creation and destruction operators 
of atoms in the radial ground state on the l-th lattice 
site. The hopping matrix element between adjacent 

sites is t = -fdzwf(z)(-^-£z+V op (zf)wi +1 (z), 

where wi(z) is a function localized on the Z-th lat- 



tice site. The on-site inter-atom interaction is U = 
g J dr± |0 o (r^)| 4 f dz\wi(z)\ 4 , where g — 4nh 2 a/m and 
a is the s-wave scattering length. 4>o(r±) is the wave 
function in the radial direction. 

Expanding around the minima of the optical lattice po- 
tential wells in a harmonic approximation, the well trap 
frequency is lo s = s 1 / 2 ^-, which is assumed to be much 
larger than the harmonic trapping frequency along the 
z-axis. Approximating the localized atomic wave func- 
tion as a gaussian at the potential minima of Z-th site, 

wi(z) = (^) 1/4 exp (-^t(z - Zl f), one obtains 
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Taking the wave function in the radial direction to be the 
ground state wave function of a harmonic oscillator, one 
obtains (d = %/k is the lattice period) 
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where a, = 



h 



We restrict ourselves to the superfluid solutions of Eq. 
lJT)l. when the phase coherence between wells is well de- 
fined in th e op tical lattice. The superfluid-Mott transi- 
tion 0, 0, only arises when the interaction energy is 
much larger than we consider here. Due to Bose conden- 
sation, the destruction operator of an atom at site I can 
be written as ai — + ifii, where = (a;) is the con- 
densate wave function and -0; is the non-condensate field 
operator. This leads to the generalized discrete Gross- 
Pitaevskii equation |l2| (we set h = 1 from now on) , 

ij^i = + +U(n c i + 2n ; )$i, (4) 

where (t) is the number of condensate atoms on the l- 
th lattice site. This includes the time-dependent Hartree- 
Fock mean field 2Uhi(t) arising from the non-condensate 
atoms on the Z-th site. Eq. (0J reduces to the usual dis- 
crete tight binding GP equation 0, @ if all the atoms 
are assumed to be in the condensate ( i.e., hi — 0). In- 
troducing phase and amplitude variables, $/ = 
Eq. |@J reduces to 



dt 



= -2t 



^n^ +1 sm(6i +1 - Bi) 
+ 2t^Jn c l n c l _ 1 sin(6» ; - 0j_i), 
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In Eq. JfJJ), e'i is an energy of a condensate atom, which 
reduces to the chemical potential of the condensate /Lt c o 
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FIG. 1: The Bogoliubov excitation energy E q (normalized to 
t) in the first Brillouin zone, for a — 1 and 6, as a function of 
the quasi-momentum q. The dashed lines give E q = cq. 

in static thermal equilibrium. The equilibrium solution 
(9i and rtf are independent of I) is 9i{t) = —[i c ot- 

The T = Bogoliubov excitation spectrum 0, 0, 
for a uniform optical lattice is easily obtained by ignoring 
the non-condensate atom term (hi = 0) and considering 

■5n?. The 



c _ „c0 . 



small fluctuations from equilibrium, n 
normal mode solutions of Eqs. 10 and JJJJ are 

69i = 59(q)e i[qld - E -> t] , 5n c i = 5n c (q)e i[qld - E " t] , 

with the Bogoliubov Bloch excitation energy, 
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At sin ^- is the kinetic energy associated with 



Here e q 

tunneling. This T = excitation spectrum is shown 
in Fig. 1 for two values of the interaction ratio a = 
Un c0 jt. For small q, the spectrum is phonon-like E q ~ 

cq, with the phonon speed c = \/2td 2 U n c0 — \J U ^»° , 

where m* = is an effective mass of atoms in the 
optical lattice. We call attention to an important feature 
of Eq. (jHJl. For a < 6, E q bends up before bending 
over, as q approaches the BZ boundary. This "anomalous 
dispersion" also occurs in superfluid 4 He |4j,[5|. For a > 
6, the spectrum simply bends over as q increases. 

As we noted earlier, the energy of a condensate atom 
el on site / is given by the r.h.s of Eq. 0. In static 
thermal equilibrium, when all sites are identical, the 
condensate chemical potential is given by /i c o = e c0 = 
-2t + U[n ca (T) + 2n°(T)]. We calculate n c0 (T) self- 
consistently by calculating the number of non-condensate 
atoms n°(T) at a site using the Bogoliubov excitations 
given by the static Popov approximation. That is, we 
ignore the dynamics of the non-condensate [i.e., we set 
Shi(t) = in Eq. ©]. This gives the Bogoliubov-Popov 
excitation spectrum E q which is identical to Eq. (JHJ, ex- 
cept that now n c0 (T) is the temperature-dependent num- 
ber of condensate atoms at any lattice site. 

Expressing n°(T) in terms of these Bogoliubov-Popov 
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FIG. 2: The site condensate fraction n c0 /n as a function of 
temperature, where s is the strength of the optical lattice 
potential. Following Ref. Q, we take n — 1200. 



excitations, we have 
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[(u 2 q + vl)f(E q ) + vl], 



(9) 



with the standard Bogoliubov transformation functions 



l/2(E q /E q + 1), 



\/2{E q /E q - 1), and 



f°(E q ) = [exp(/3E q ) - I]" 1 . The Hartree-Fock (HF) ex- 
citation spectrum is E q = [e q + 2U(n c0 + h )] — £iq — 
4t sm 2 (^)+Un c0 . Apart from the limiting case of a <C 1, 
we note that E q is quite different from the HF spectrum 
E q . Strictly speaking, in our ID model, there are no so- 
lutions of Eq. @ for an infinite system because the sum- 
mation diverges due the contribution from small momen- 
tum, in accordance with the Mcrmin- Wagner- Hohcnbcrg 
theorem. However, we consider finite systems which in- 
troduces a lower momentum cutoff q c = 2ir/Id. (As a 
typical value , we take the total number of lattice sites 
/ = 250.) In this case, Eq. @ has a solution and the 
condensate can exist. We are only considering a ID op- 
tical lattice for illustration and emphasize that our main 
conclusions also hold for 2D and 3D optical lattices 

The condensate number n c0 (T) is found by solving Eq. 
© self-consistently for a fixed value of n. In Fig. 2, we 
plot n c0 (T)/n as a function of T. Using these results, 
we plot the dimensionless interaction a = Un cQ (T)/t as 
a function of temperature in Fig. 3. The ratio a and 
the results in Fig. 3 will play a crucial role in our subse- 
quent analysis. Since we limit our discussion to the first 
energy band of the optical lattice, our results only apply 
when s ^> 2k^TjE]i. Higher bands would be thermally 
populated if we considered lower values of s. 

We next turn to the calculation of Landau damping. 
A collective mode in a condensate coupled to a thermal 
cloud of excitations in thermal equilibrium has a complex 
frequencVj u> = E q — \T q , with the Landau damping given 
by HQ 

T q = it Y, \M q . M \ 2 [f (E p ) - f (E q+p )] 

H>*?c 

x 5 (E q + E p - E q+P ) . (10) 



FIG. 3: Plot of the interaction parameter a = Un c0 (T) /t, as 
a function of temperature, for several values of s. The number 
of condensate atoms at a lattice site n c0 (T) is given in Fig. 
2. The dashed line is a — 6. 

Here E p is the Bloch Bogoliubov excitation energy given 
in Eq. (|5J) for a uniform condensate in an optical lat- 
tice, and the momentum sum is over the first Brillouin 
zone. The expression in Eq. (jl0|) describes a condensate 
excitation of (quasi) momentum q being absorbed by an 
excitation p of the optical lattice thermal gas, leading to 
a thermal excitation with momentum q + p. 

The energy conservation condition E q + E p — E q+p 
needs to be satisfied in Eq. (|10fl . This is illustrated in 
Fig. 4 16]. Clearly, the intersection at {q + p, E q+p ) re- 
quires that the dispersion relation E q first bends up as q 
increases, before bending over. This anomalous disper- 
sion is also a feature of phonons in superfluid 4 He |4j, |5j . 
The values of (q, p) satisfying E q + E p = E q+p is shown 
in Fig. 5. For a given q, we see that as a — > 6, the value 
of p decreases to zero. There is no solution for a > 6, 
indicating the disappearance of Landau damping. 

The matrix element in Eq. I|10(l is given by [l4l fl5| 



M q pi - P2 2U y j / t [(^pi^p2 ~^ ^pi^pi v pi u p2 ^)u q 
G„ 

-(u pl u P2 +v Pl v P2 - u Pl v P2 )v q ]S p+Pu P2+G „, (11) 

where G n = is a reciprocal lattice vector (n is an 
integer). In Fig. 5, the curves of the solution of the 
energy conservation condition never go across the dashed 
line q + p = ir/d. Thus Umklapp scattering processes 
(G n 7^ 0) do not contribute to Landau damping [lfij . 

When the excitation q has a wavelength much larger 
than the thermal excitation p (i.e. q <C p, ir/d), the en- 
ergy conservation condition E q + E p = E q+p in Eq. I|10|) 
means that for q <^ p, the Landau damping of the excita- 
tion E q = cq comes from absorbing a thermal excitation 
E p with a group velocity equal to c ^^] . The unique 



value of po such that 



dE p 
dpo 



c is found to be given by 



the condition 4 sin 2 (\p d) = -(a - 2) + y/2(a + 2). This 
is only valid for a smaller than 6, such that po q. 
For q <^ p, the matrix element reduces to (the analogue 
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FIG. 4: The Bogoliubov excitation energy E q for a < 6. The 
intersection of the two dispersion curves at p + q shows that 
the energy conservation condition is satisfied. 




Beliaev damping (present even at T — 0) involving the 
spontaneous decay of an excitation into two excitations. 
All such processes also involve an energy conservation 
factor of the kind 8(E\ — E% — E3), which can only be 
satisfied if the excitations exhibit anomalous dispersion, 
i.e., a < 6D, where D is the dimension of the optical 
lattice. Finally, for a fixed value of a < 6D, there is a fi- 
nite threshold momentum q* such that the decay (Beliaev 
damping) of an excitation E q is impossible when q > q*. 
The same phenomenon has been discussed for phonons 
in superfluid 4 He 0, 0, 0] and one finds that q* ~ p 
Q . This disappearance of Beliaev damping for q > q* at 
T = should be easy to confirm experimentally. 

In conclusion, assuming that only the lowest energy ex- 
citation band is thermally occupied, we have shown that 
Landau damping of Bogoliubov excitations in uniform 
optical lattices is only possible if they exhibit anomalous 
dispersion. 
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FIG. 5: The values (q,p) satisfying the energy conservation 
condition E q + E p = E q + P , for several values of a. 



of Eq. (10.76) in Ref. U\ ) 
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where iV c0 = n c0 I. In deriving Eq. ifT^ from Eq. (fTTfl ■ 
we have used the energy conservation factor in Eq. (|10|l . 
Using the delta function for energy conservation to inte- 
grate over p, the Landau damping is given by 
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One finds there is no damping when a > 6 (see Fig. 3 
for a as a function of s and T). T q diverges at a — 6, but 
this is due to the ID nature of our system and does not 
occur in 2D lattices. Fig. 3 shows that, for n = 1200 0, 
a is smaller than 6 only when T is very close to T c , for s 
in the range 6 ~ 9. For a smaller number of atoms on a 
lattice site, the temperature range where a is less than 6 
becomes larger. 

For our ID model to apply, one would need a much 
tighter magnetic trap (in the radial direction) than used 
in Ref.pJ. A 3D optical lattice (without any magnetic 
trap) might allow the clearest check of our predictions of 
the disappearance of damping. 

One also has damping from collisions which transfer 
atoms between the condensate and thermal cloud Il7l and 
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